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Abstract. We show that by changing a single non-dimensional number, the thermal 
Rossby number, global atmospheric simulations with only axisymmetric forcing pass 
from an Earth-like atmosphere to a superrotating atmosphere that more resembles the 
atmospheres of Venus or Titan. The transition to superrotation occurs under conditions 
in which equatorward-propagating Rossby waves generated by baroclinic instability at 
intermediate and high latitudes are suppressed, which will occur when the deformation 
radius exceeds the planetary radius. At large thermal Rossby numbers following an ini- 
tial, nearly axisymmetric phase, a global baroclinic wave of zonal wavenumber one gen- 
erated by mixed barotropic-baroclinic instability dominates the eddy flux of zonal mo- 
mentum. The global wave converges eastward zonal momentum to the equator and de- 
posits westward momentum at intermediate latitudes during spinup and before super- 
rotation emerges, and the baroclinic instability ceases once superrotation is established. 
A global barotropic mode of zonal wavenumber one generated by a mix of high- and 
low-latitude barotropic instability is responsible for maintaining superrotation in the 
statistically steady state. At intermediate thermal Rossby numbers, momentum flux 
by the global baroclinic mode is subdominant relative to smaller baroclinic modes, and 
thus strong superrotation docs not develop. 



1. Introduction 

Superrotation 1 is a feature of a number of planetary at- 
mospheres, although the mechanisms that give rise to it may 
differ from case to case. It is a well-known result that an 
axisymmetric atmosphere cannot superrotate if there is a 
diffusive mechanism that acts to mix angular momentum 
down-gradient [Hide, 1969]. Superrotation can then only 
arise if there are non-axisymmetric eddy motions that act to 
transfer momentum upgradient and into an angular momen- 
tum maximum. Rossby waves are a potential mechanism for 
this: it is a property of such waves that momentum is trans- 
ported in the opposite direction to their group velocity, and 
thus momentum converges in their source region [Thompson, 
1971]. This mechanism gives rise to the eastward surface 
winds in midlatitudes in the Earth's atmosphere (where 
baroclinic instability is the source) and may lead to superro- 
tation in two-layer models of Earth's atmosphere provided 
the tropical forcing is sufficiently non-axisymmetric [Suarez 
and Duffy, 1992; Saravanan, 1993]. Equatorial superrotation 
is present on Jupiter and Saturn [Flasar, 1986], requiring 
some form of up-gradient angular momentum transfer [Read, 
1986, e.g.]. Rossby waves may also be responsible for equato- 
rial superrotation on the giant planets (where convection is 
likely to be the source) [Showman, 2007; Scott and Polvani, 
2008; Schneider and Liu, 2009; Lian and Showman, 2009], al- 
though another mechanism involving equatorial Kelvin waves 
has been proposed [Yamazaki et al., 2005]. Zonal jets and 
equatorial superrotation on Jupiter and Saturn may alter- 
natively result from mechanisms involving deep convection 
[Heimpel and Aurnou, 2007; Kaspi et al., 2009] which are 
unlikely to occur in shallow, terrestrial atmospheres. 
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Without considering the origin or nature of the eddy mo- 
tion, Gierasch [1975] posited that large-scale eddies could be 
parameterized by a term in the zonal momentum equation 
proportional to the angular velocity gradient, and so with the 
same form as a molecular viscosity. The effect was further 
explored by Rossow and Williams [1979] and others in the 
context of superrotation on Venus. Such an eddy viscosity 
would tend to bring the fluid into solid body rotation, po- 
tentially leading to superrotation when coupled to a steady 
overturning circulation. This effect has been invoked, at 
least in part, by a number of subsequent studies as being 
the cause of superrotation on Venus and Titan [Rossow and 
Williams, 1979; Del Genio and Zhou, 1996; Hourdin et al., 
1995; Luz et al, 2003; Lee et al, 2007], with the proposed 
source of the eddies being off-equator barotropic instability. 
Williams [2003] discussed how the eddy fluxes associated 
with barotropic instability can give rise to superrotation, 
and noted that tropical instabilities may also play a role. 
Yamamoto and Takahashi [2004; 2006] found a similar effect, 
with the eddies in their simulations being at least partially 
maintained by thermal tides. 

In this study our focus is on the underlying dynamical 
mechanisms producing superrotation and its parametric rela- 
tionship to Earth's atmospheric circulation. Thus, we use a 
three-dimensional model with a fairly simple treatment of the 
diabatic physical processes, but explore a dynamically wide 
parameter regime with a goal of passing from an Earth-like 
atmosphere to a superrotating one with a change in a single 
nondimensional parameter, so exposing the mechanism as 
clearly as possible. In section 2 we describe the governing 
equations and numerical model. In section 3 we identify the 
key nondimensional numbers, and in section 4 we describe 
numerical experiments. We follow this, in section 5 by a dis- 
cussion of the dynamical processes producing superrotation 
and, in section 6, of the relevance to planetary bodies in the 
Solar System. We conclude in section 7. 
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2. Numerical Model Description 

Our numerical model integrates the dry primitive equa- 
tions of motion of an ideal gas on a sphere using a spectral 
method [Gordon and Stern, 1982] at T42 resolution (test 
cases were performed at higher and lower resolutions with 
few essential differences found). The main differences from a 
comprehensive GCM are that our model has no condensate 
(for example water vapor for Earth, methane for Titan) and 
that the radiative effects are represented by a Newtonian 
cooling. These simplifications are severe, but allow us to 
focus on the dynamical aspects of the circulations of the 
systems. The simplest formulation of diabatic atmospheric 
processes is to use linear relaxation, or Newtonian cooling, 
toward a specified temperature profile. The forcing profile 
is fixed and a timescale, the radiative cooling time, param- 
eterizes the time dependence. Our approach follows Held 
and Suarez [1994] with some differences in detail, as we now 
describe. 

The forcing profile is assumed to be zonally symmetric, 
and it is assigned a fixed latitude dependence with a non- 
dimensional parameter, Ah, specifying the cquator-to-pole 
temperature gradient, T a = T[l + Ah/3(1 — 3 sin ip)] with 
lowest-level temperature T and global average surface tem- 
perature T. We obtain the vertical structure of the forcing 
by assuming a relaxation to vertical temperature profiles 
given by a moist adiabat with a lapse rate of 6 K/km, thus 
our radiative-equilibrium profile is generally stable to dry 
convection. The forcing profile is capped with an isothermal 
stratosphere by not allowing temperatures to drop below 
70% of T . The structure of the forcing profile is shown in 
Fig. 1 (dotted contours). The specification of our radiative 
forcing is completed by setting the radiative relaxation time, 
and we use values from 40 days in the free troposphere to 
4 days in the boundary layer, the top of which is fixed at 
p/Ps = 0.7. 

Frictional stresses exchange angular momentum between 
the surface and atmosphere, and we specify this mechanism 
as a linear Rayleigh friction in the lower atmosphere. The 
Rayleigh relaxation time is taken to be 1 day at the surface 
and increases to the top of the boundary layer, above which 
there is no explicit friction except for small amount of vertical 
diffusion, with a diffusion coefficient of v = 0.01 m 2 /s, that 
acts to smooth gridscale noise. A fourth-order hyperviscosity 
is applied to dissipate energy at the grid scale. 

3. Nondimensionalization and Parameter 
Regimes 

Given the relatively simple thermodynamic and frictional 
schemes described above it becomes possible to determine 
the nondimensional parameters that govern the behavior of 
the system. We begin with the primitive equations subject 
to Rayleigh friction and Newtonian cooling. The horizontal 
momentum equations may be written in vector form [Valhs, 
2006, Section 2.2.5] as 



du 

~dt 



+ v-\7u + fxu = -Vp$ 



(3.1) 



where v = Mi + uj + wk, u = ui + vj and / = /k = 20smp]i., 
with (i,j,k) being the unit vectors in the zonal, meridional 



Table 1. Parameters for experiment design. A dash indicates 
the same value as that of the entry on its left. 



Rot 



0.02 



1.3 



10.5 



~o 6.4 x 10 e m 

n 7xl0- 5 s~ 1 

T 285 K 

A H 0.2 



8 x 10 5 m 2.8 x 10 5 m 



and local vertical directions. The hydrostatic, thermody- 
namic and mass continuity equations complete the model. 

Denoting non-dimensional variables with a hat, we non- 
dimensionalize the above equations by setting (u, v) = 
(u,v)/U; w = wa/(UH); $ = &/(2f2Ua); t = 2(2t; f = 
f/(2Q) = simp, and (0,f) = (0, T)/(T A H ); with U a 
characteristic horizontal wind speed, a the planetary ra- 
dius, H the scale height and i? the rotation rate. For the 
scaling of winds we use the thermal wind relation, whence 
U — RT Ah /(2fla). Making these substitutions into the 
above dimensional primitive equations we obtain 



du 



+ Ro T (v ■ X7u) + / x u = -Vp$ - Eu . (3.2) 



where the two nondimensional parameters are 

r 



An Ekman number: E 



2(2' 



(3.3a) 



A thermal Rossby number: Rot 



U _ (RToA H ) 



2f2a (2f2a) 



(3.3b) 



A third non-dimensional parameters appears in the ther- 
modynamic equation, namely the nondimensional radiative 
relaxation time, f = 2J?t. Other nondimensional parameters 
of the model not explicit in this derivation which we do not 
vary include the Prandtl number and the Rayleigh number. 
The aspect ratio, H/a with scale height H, enters implicitly 
through the fixed depth of the boundary layer. Note that 
geostrophic balance, which is implied in the thermal wind, 
may not be a good approximation at high Rossby number 
so that although Rot remains a valid non-dimensional pa- 
rameter it may not be an accurate measure of the Rossby 
number itself. The thermal Rossby number has been used 
as a control parameter in rotating annulus experiments [e.g., 
Hide, 1958; Geisler et al, 1983] as well as numerical studies 
of Earth's atmosphere [e.g., Held and Hou, 1980]. 

In order to focus on the emergence of superrotation, we 
report only on the effects of the thermal Rossby number. The 
Ekman number determines the strength of the drag at the 
bottom, and so the surface wind, and is not responsible for 
the qualitative difference in the regimes. The radiative relax- 
ation time in part determines the strength of the Hadley cell 
[Held and Hou, 1980], but we will not vary this parameter. 

4. Numerical Experiments 

We now describe the results of a set of numerical experi- 
ments aimed at characterizing the transition of our model 
atmosphere from an Earth-like case to one which superro- 
tates. By definition, superrotation occurs when the axial 
angular momentum of the atmosphere exceeds solid-body 
rotation, and this could in principle occur at mid-latitudes 
without prograde equatorial winds. However, this would 
require angular momentum of the atmosphere to increase 
poleward, and it is likely this situation would be unstable to 
symmetric perturbations. An adjustment process by sym- 
metric eddies thus would eliminate the off-equatorial angular 
momentum maximum. Superrotation is virtually always 
characterized by prograde equatorial winds and for all our 
experiments prograde equatorial winds, if present, mark the 
atmospheric angular momentum maximum. 

It is customary to describe Titan and Venus as being in a 
'slowly rotating' regime compared to Earth and, for Venus, 
the slow rotation compared to Earth is the main factor giving 
it a large thermal Rossby number. Titan, however, is also 
substantially smaller than Earth, and this also contributes 
to its large thermal Rossby number. In the present study we 
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Figure 1. Thermal forcing profile used for the experiments (dotted) and the zonal and time mean 
temperatures over the last 360 days of 1080 day integrations (gray line) for each of our three model 
cases. Contours are spaced evenly at intervals of Ah/4 from 0.7T o to T , where T is the global mean 
temperature. 



Ro = 0.02 Ro = 1.3 Ro = 10.5 




latitude latitude latitude 

Figure 2. Zonal and time mean mass streamfunction (shaded) and zonal winds (contours, zero-line bold) 
for Rot = 0.02, 1.3, and 10.5. The top row displays the full three dimensional simulations and the bottom 
row displays axisymmetric simulations of the same suite of parameters. Zonal wind contours begin at 
m/s (bold line) and are spaced by 10 m/s increments. Mass streamfunction contour ranges are reduced 
by a factor of 10 for each increase in Ro, and the maximum mass flux for each case is marked in white 
text [Tg/s]. 



vary only the thermal Rossby number. In practice, because 
the numerical code is dimensional, we must vary a dimen- 
sional parameter and we choose to vary only the planetary 
radius and not the rotation rate. By making this choice 
the Ekman number and the nondimensional radiative relax- 
ation time remain fixed. Both of these parameters are in 
fact quite different for Titan, but the simplification is useful 
in providing clean numerical experiments. Table 1 outlines 
the dimensional parameters in our experiments, focusing on 
three representative cases, Rot = 0.02, 1.3, and 10.5. 

4.1. The standard cases 

Global, zonal-mean diagnostics from the three standard 
simulations are shown in Fig. 1 and Fig. 2. Fig. 1 displays 
the zonal and time mean temperature structure (gray con- 



tours) and the forcing profile (dotted contours) averaged over 
the final 360 days of 1080 day integrations. Contours are 
spaced evenly at intervals of Ah/4 from 0.7T o to T with 
T the global mean surface temperature. The most marked 
difference between the three cases is that as Rot increases, 
the latitudinal extent of horizontally uniform temperature 
increases. This feature marks the expansion of the tropical 
circulation pattern; in the Rot = 10.5 case, there is only a 
small remnant of the original horizontal temperature gradi- 
ent, now isolated to regions poleward of 60° N/S latitude. 
Thus the state of the atmosphere becomes more barotropic 
as Rot increases. 

Fig. 2 shows the time and zonal mean overturning circula- 
tion (shaded) and zonal winds (contours) for our numerical 
experiments. The top row displays the full three-dimensional 
simulations while the bottom row shows axisymmetric ver- 
sions of the same three cases. The latter were obtained by 
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initializing the three-dimensional model without initial seed 
perturbations, and were verified to be very nearly axisym- 
metric throughout the simulations. The color scale of the 
overturning circulation has been reduced by a factor of 10 
between successive increases in Rot while the scale of zonal 
winds is fixed (10 m/s contour spacing; zero- wind line bold). 
Between Rot = 0.02 and Rot = 1.3, the Hadley cell for 
the three-dimensional simulations widens and substantially 
weakens while the indirect cell due to eddy motions in the 
mid latitudes is pushed poleward. The Hadley cell in the 
Rot = 10.5 case is substantially weaker aside from a vis- 
cous cell contained in the boundary layer (pressures greater 
than p/ps = 0.7). The mean circulation above the boundary 
layer in this case slants upward and poleward away from the 
equator. This slant structure is characteristic of symmet- 
ric instability arising from the misalignment of contours of 
constant angular momentum and isentropes. 

The zonal mean zonal winds in the top row of Fig. 2 
transition from an Earth-like, non-superrotating regime at 
Rot = 0.02 to a fully supcrrotating regime at Rot = 10.5. 
(Surface zonal winds in the Rot = 10.5 case are all easterly 
because of the use of a Rayleigh drag in the lowest few model 
layers and not just at the ground; however, the integrated 
torque over the boundary layer is zero, as is required for a 
steady-state.) The transition between regimes occurs near 
the Rot = 1.3 case, in which the zonal wind has Earth-like 
isolated mid-latitude jets and also equatorial superrotation. 
This is to be expected because, as Ro increases past unity, 
the scale of baroclinic turbulence (which has a significant 



impact on the zonal wind structure) can no longer fit within 
the confines of the spherical geometry of the planet as also 
noted by Williams [2003]. 

To see what causes baroclinic instabilities to become inef- 
ficient, compare the deformation radius (which largely deter- 
mines the scale at which baroclinic cyclogenesis occurs) to 
the domain size. The deformation radius, Ld = NH/f, does 
not significantly deviate from that of our thermal forcing 
profile in any of our simulations. Therefore Ld is essen- 
tially fixed; with typical mid-latitude values of H ~ 10km, 
/ ~ n, and iV 2 ~ g/T (g/Cp + dT a Jdz) ~ lO^s" 2 we 
find Ld ~ 1500 km ~ aEarth/4. Also note that, from (3.3b), 
the thermal Rossby number Rot is proportional to a~ 2 , and 
since in varying Rot we are changing the radius of our planet 
(while keeping the deformation radius fixed), we expect baro- 
clinic instability to be severely limited for cases with a < Ld, 
or in other words for Ro > 16-Rot, Earth = 0.32. Figure 3 
shows that in the statistically steady state in the upper tro- 
posphere, the local Rossby number, which we estimate as Ro 
fa u/ (| f\ a cos ip) (that which is realized by the simulation), 
only moderately and locally exceeds unity in the Rot = 1.3 
case, indicating that baroclinic instability is still able to oc- 
cur. In the Rot = 10.5 case the local Rossby number reaches 
values of Ro = 2-5, and in this case baroclinic instabilities 
are inhibited because the deformation radius, and hence the 
scale of the instability, is becoming comparable to or larger 
than the planetary radius. 




Figure 3. An estimate of the local Rossby number, Ro « w/(|/|acos ip) with zonal mean zonal wind u, 
(absolute value of the) Coriolis parameter |/|, planetary radius a, and latitude ip for our three standard 
cases averaged over the last 360 days of 1080 day integrations. Shaded contours are on the same color 
scale in each panel, and are spaced from to 3 by 0.5. The line spacing changes for each panel: Rot = 
0.02 lines are spaced from to 0.1 by 0.02, Rot = 1.3 are spaced from to 2 by 0.5, and Rot = 10.5 are 
spaced from to 5 by 0.5. The bold line marks Ro — 0. 



The relative weakness of baroclinic turbulence does 
not guarantee the development of superrotation, since 
there must also be a source of momentum convergence 
at the equator. However, the presence of baroclinic 
instabilities does not of itself preclude superrotation ei- 



ther, since in the absence of an equatorial critical layer, 
Rossby waves generated in midlatitudes pass through 
the tropics without breaking and so without decelerat- 
ing the flow (see schematic in Fig. 8). The mechanism 
leading to superrotation at the equator is the topic of 
the Section 5. 
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Figure 4. Zonal and time mean zonal winds (solid line, m/s) and potential vorticity (dotted line) for the 
steady-state of Rot = 0.02, 1.3, and 10.5 and at the 200 (top) and 500 (bottom) hPa levels. 



Fig. 4 displays the zonal- and time-mean zonal winds 
(solid lines) and potential vorticity (PV, dotted lines). 
The panels are arranged similarly to Fig. 2, with Rot 
increasing from left-to-right, while the top and bot- 
tom rows display these diagnostics at 200 hPa and 
500 hPa, respectively. Especially in the middle atmo- 
sphere, the Earth-like case with Rot — 0.02 has very 
non- uniform PV, while at Rot = 1.3 the PV becomes 
nearly piece- wise homogeneous, forming staircase-like 
structures. The Rot — 10.5 case has very non-uniform 
PV at all levels. 

Fig. 5 displays a snapshot on day 1080 of geopotcntial 
height on the 400 hPa level in our three standard cases. 
Contours are spaced evenly from 900 to 980 m in incre- 
ments of 10 m. The Rot =0.02 case is evidently domi- 
nated by disturbances with zonal wavenumbers ^4 or 
larger. These disturbances exhibit a northeast-southwest 
tilt in the northern hemisphere (northwest-southeast in 
the southern hemisphere) as is characteristic for baro- 
clinic Rossby waves. The Rot =1-3 and Rot =10.5 show 
clear evidence for a zonal wavenumber one disturbance 
with coherence at all latitudes. 

4.2. The axisymmetric cases 

The axisymmetric simulations in the bottom row of 
Fig. 2 do not develop superrotation (as they cannot in 
an axisymmetric flow), aside from very weak equatorial 
winds at high Ro. (The very weak superrotation is due 
to a small amount of numerical mixing.) It is clear 
from these plots that strong superrotation requires the 
influence of non-axisymmetric eddies. Subtropical jets 
form due to the advection of high-angular-momentum 
air poleward from the equator, and the latitude of the 



jets increases with increasing Rot as expected from 
axisymmetric, inviscid theory [Held and Hou, 1980]. 

As Rot is increased, the axisymmetric case develops 
a wider swath of time- dependent symmetric overturning 
that, in the time-mean, has the appearance of a Hadley 
circulation (especially in the Rot = 10.5 case). This 
axisymmetric component does not have a significant 
projection on the three-dimensional simulations for the 
Rot = 1.3 case. For Rot — 10.5, the axisymmetric 
circulation is clearly visible in the time-mean Hadley 
circulation. The similarity in the magnitude of the ax- 
isymmetric overturning in the Rot = 10.5 case relative 
to the full three-dimensional simulation suggests axisym- 
metric theory for Hadley cells is a better approximation 
at higher Rot- 

4.3. The sponge cases 

Previous studies of superrotation have often at- 
tributed equatorial eddy momentum convergence to 
high-latitude barotropic instabilities feeding off of hor- 
izontal shear of the jets (see e.g. Williams [2003] or 
Hourdin et al. [1995]). However, the eddies that produce 
the superrotation cannot be Rossby waves emanating 
from high-latitudes; if such waves were to radiate from 
a high-latitude source region and break at low latitudes, 
they would deposit westward momentum thus decelerat- 
ing equatorial winds. We now demonstrate in a set of 
test simulations that the low-latitudes play an equally 
important role in genesis of the disturbances that give 
rise to superrotation. 

In these test simulations, we impose a damping, or 
"sponge" on the non-axisymmetric component of all fields, 
but only in a select band of latitudes. The sponge takes 
the form of a linear damping of the non-axisymmetric 
component at a timescale of 10 minutes, which results 
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Figure 5. Snapshots at day 1080 of the 400 hPa-level geopotential height in our three standard cases. 
Contours are spaced evenly from 900 to 980 m in increments of 10 m. 




Figure 6. Hadley cell and zonal mean zonal winds for test simulations of the Rot = 10.5 case with 
sponge layers (contouring is as described in Fig. 2 for the Rot = 10.5 case). The left panel displays the 
control case with no sponge applied. A sponge has been applied to latitudes poleward of 45 N/S latitude 
in the middle panel. The right panel has a sponge applied at the closest two grid cells to the equator, 
corresponding to the area between 3 N/S latitude. 



in an essentially axisymmetric state in the latitude band 
where the sponge is applied. These test simulations 
are aimed at understanding the influence of equatorial 
vs. high-latitude eddies on the zonal mean zonal flow. 
The first of these cases labeled "Hi-lat sponge" in Fig. 
6 has the sponge applied to regions poleward of 45N/S 
latitude. Zonal mean zonal winds in this high-latitude 
sponge test case are only weakly superrotating at the 
equator, lending evidence of a high-latitude influence on 
the development of equatorial superrotation. The second 
case labeled "Lo-lat sponge" in Fig. 6 has the sponge ap- 
plied to one grid-cell on either side of the equator, which 
at T42 resolution corresponds to roughly 3N/S latitude. 
It is remarkable that in this low-latitude sponge case, the 
flow is nearly identical to the axisymmetric simulation 
at all latitudes (compare with lower, right-hand panel of 
Fig. 2), establishing that the low-latitudes play a crucial 
role in the development of superrotation in the Rot = 
10.5 case. 



5. Dynamical Interpretation 

We now apply several diagnostics to our simulations 
and give an interpretation of them. We begin by com- 
paring the Eliassen-Palm fluxes between the Ro^ — 
0.02 and Rot — 10.5 cases. We then use a space-time 
fourier analysis of the Rot = 10.5 case to identify the 
disturbances that give rise to superrotation. We end 
with a characterization of the disturbances contributing 
to the zonal mean zonal momentum and describe the 
mechanism by which superrotation is established and 
maintained. 

5.1. Eliassen-Palm fluxes 

In Fig. 7, we show the (full ageostrophic) Eliassen- 
Palm (EP) flux (arrows) and its divergence (shaded) 
during representative times of the Rot = 0.02 and 10.5 
simulations. The EP flux is often interpreted as a mea- 
sure of the flux of wave activity, and its divergence gives 
the acceleration of zonal momentum due to eddy pro- 
cesses [Edmon et at, 1980; Vallis, 2006]. In the figures, 
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Figure 7. EP fluxes and their divergences. The left panel displays the Rot = 0.02 (Earth-like) case after 
it has reached a statistically steady state (during days 1000 to 1080). The middle panel displays the Rot 
— 10.5 case (Titan-like), during active spinup, of the layer near p/p s = 0.4 (days 215 to 305), and the 
right panel displays the Rot = 10.5 case after reaching statistically steady state (days 1000 to 1080). 



Earth-like Superrotating 




Figure 8. Schematic representation of eddy-mean-flow interaction in the Earth-like (left) and Titan-like 
(right) regimes. "Wave source" regions are intended to represent the location in latitude and characteristic 
phase speed of disturbances relative to the zonal wind. In so far as dissipative effects are otherwise small, 
meridional propagation of these disturbances continues until either a critical layer is reached, denoted 
by arrows intersecting the zonal wind, or they are reflected by the vorticity gradient, denoted by arrows 
bounded by horizontal lines. (The generation of Rossby waves at the equator in the Earth's current 
atmosphere is in fact relatively weak.) 



dark (light) colors indicate acceleration (deceleration) 
of the zonal winds. The left panel in Fig. 7 shows 
the EP flux and its divergence for the Earth-like case, 
Rot = 0.02, in the statistical equilibrium phase. This 
has the characteristic pattern of mid-latitude baroclinic 
turbulence; the turbulence accelerates surface winds in 
mid-latitudes, decelerates the upper-level winds, and 
radiates Rossby waves equatorward where they are ab- 
sorbed and decelerate subtropical winds. The middle 



panel displays the acceleration phase of the equatorial 
winds of the Rot = 10.5 case, during which time there 
is a divergence of EP fluxes at the 400 and 700 hPa 
levels (p/ps = 0.4 and 0.7). There is no evidence for a 
divergence of vertical EP fluxes. Instead, a disturbance 
is contributing acceleration of equatorial zonal winds 
and deceleration of mid- and high-latitude zonal winds 
through a horizontal flux of momentum. 

Once the Rot — 10.5 simulation is spun-up very little 
EP flux activity occurs in the upper, superrotating levels 
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as seen in the right panel of Fig. 7. There is significant 
eddy activity at higher latitudes but this clearly does not 
directly affect the equatorial regions aloft. The relative 
lack of eddy activity or mean flow (see top right panel of 
Fig. 2) in the upper, superrotating levels of the equator 
is similar to the finding of Saravanan [1993] that once 
equatorial superrotation is established the absence of 
equatorial torques allows the superrotation to be main- 
tained for very long periods with little equatorial eddy 
activity to provide a momentum convergence. After the 
spinup phase, there need be only very occasional eddy 
momentum fluxes to sustain superrotation, and this is 
the case in our Rot = 10.5 simulation. 

To summarize, (and see also Fig. 8) in the Earth-like 
case, Rossby waves, generated in mid-latitudes by baro- 
clinic instability, propagate upwards and equatorwards 
and break, decelerating the tropical flow and inhibiting 
any superrotation that might arise from other mecha- 
nisms. In the superrotating case, Rossby waves from 
mid- or high-latitude regions either do not exist or do 
not interact with the mean flow at low latitudes. 

5.2. Eddy momentum flux cospectra 

We now perform spectral diagnostics of wave activity, 
primarily of the Rot — 10.5 case, in order to establish 
those disturbances which contribute to the acceleration 
of the mean flow. Because we are interested in the in- 
fluence of eddy momentum convergence on zonal mean 
zonal winds we compute cospectra of the horizontal flux 
of zonal momentum given by 

K niU = 2<Re(L/; ;( X:J> ( 5J ) 

Here, (If, V) are the fourier transforms of the non- 
axisymmetric components of (u,v), starred quantities 
are complex conjugates, "Re" refers to the real compo- 
nent, and angle brackets denote time averaging over 
a given window. Fourier transforms are performed 
by direct FFT in space and time, resulting in fields 
gridded in zonal wavenumber, n, and frequency, ui 
[Hayashi, 1971]. Positive (negative) frequencies rep- 
resent westward- (eastward-) propagating waves. Time 
averaging is accomplished in Fourier space by smoothing 
frequencies with a gaussian window (see Randel and 
Held [1991]). 

5.2.1. The dominance of zonal wavenumber 1 

In order to identify which wavenumbers interact with 
the zonal mean zonal wind, the left panel of Fig. 9 dis- 
plays the horizontal divergence of the eddy flux cospec- 
trum dK n ^/dy at the 400 hPa level and at the equator 
for the entire first model year of the Rot = 10.5 case 
(negative divergence indicates eastward acceleration of 
the mean flow). It is clear that zonal wavenumber one 
is dominating the divergence of eddy momentum fluxes, 
primarily at frequencies of -0.1 days -1 and -0.6 days -1 . 

The middle panel of Fig. 9 displays the divergence 
of K niU at the 400 hPa level during days 250 to 260 of 
the Rot = 10.5 case, which is the most active phase 



of spinup. Wavenumber one disturbances again clearly 
dominate the divergence of eddy momentum flux, pre- 
dominantly at frequencies of —0.25 days -1 and —0.6 
days -1 . The right panel of Fig. 9 shows the divergence 
of K n>u at the 400 hPa level after the model has spun 
up, days 1080-1200. A wavenumber one disturbance is 
dominant at a frequency of -0.25 days -1 . Eddy momen- 
tum flux is concentrated in disturbances with negative 
frequencies, implying eastward propagation. However 
a substantial Doppler shift is present due to the strong 
(and latitude-dependent) zonal mean zonal wind. We 
must map frequencies of the disturbances to phase speeds 
to understand how their phase speeds compare with the 
mean flow speed. 

5.2.2. Wave-mean-flow interaction 

Linear theory predicts non-axisymmetric disturbances 
break and deposit their pseudomomentum into the zonal 
mean flow where their zonal phase speed, c, is roughly 
equal to the background zonal wind, u; the breaking 
region is referred to as a critical layer. For meridion- 
ally propagating disturbances, critical layers are ori- 
ented along latitude circles where the zonal mean zonal 
wind u w c. It is therefore instructive to derive the 
wave activity as a function of phase speed and lati- 
tude by converting the eddy momentum flux cospectra 
from wavenumber-frequency space to wavenumber-phase 
speed space and to sum the activity over wavenumber 
at each latitude. The result can then be directly com- 
pared with the zonal mean zonal wind and inspected for 
critical layers. 

The maintenance of the mean zonal winds in the 
Earth-like case by eddy-mean-flow interaction is reason- 
ably well understood and is illustrated schematically, 
along with a superrotating case in Fig. 8. In the Earth- 
like regime (Rot = 0.02), baroclinic turbulence originat- 
ing in the mid-latitudes deposits westward momentum 
in the equatorial region, decelerating the zonal winds 
there. Equatorial waves are trapped in the waveguide 
formed by the large-scale vorticity gradient, and so they 
cannot redistribute angular momentum through hori- 
zontal propagation [Gill, 1982]. If strong equatorial 
superrotation is present, the equatorial region becomes 
largely transparent to disturbances that are propagating 
equatorwards from higher latitudes, because the strong 
eastward zonal winds prevent there being a low-latitude 
critical layer, and so there is no eddy-damping of the 
eastward flow. However, and as suggested by Fig. 7, to 
generate the superrotation in the first instance there 
must be a source of wave activity at or near the equator, 
and we now look at its spectral properties. 

Fig. 10 displays the eddy momentum flux cospectra 
after converting frequency to phase speed, K n ^ Cp (inter- 
polating onto a regular grid in phase speed, as in Randel 
and Held, 1991) at the 400 hPa level and in 20-day win- 
dows centered on the day indicated (shaded, with dark 
or light indicating positive or negative values) with the 
average zonal mean zonal wind during this time overplot- 
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Figure 9. Divergence of the eddy momentum flux cospectrum, dK n ^jdy at the 400 hPa level and at 
the equator for Rot = 10.5 averaged over the days indicated. Contour ranges are the same in each panel 
(their magnitudes are arbitrary). 
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Figure 10. Phase speed-latitude plot of eddy momentum flux cospectra at the 400 hPa level (shaded) 
and zonal mean zonal wind (line) for the Rot = 10.5 case over a sequence of 20-day intervals during an 
active phase of spinup. Dark (light) regions correspond to positive (negative) eddy momentum flux. 



ted (line) at a sequence of times during active spinup 
of the Rot = 10.5 case. In the first snapshot - day 
230 - there are two global waves present with oppositely 
oriented momentum fluxes at low-latitudes. The slower 
wave with an equatorial phase speed of a ~ 2 m/s (cor- 



responding to a frequency of -0.1 days -1 ) is a westward- 
propagating Rossby wave, i.e., it is geostrophic (not 
shown), initiated at high-latitudes. This disturbance 
decelerates equatorial winds by depositing westward (ret- 
rograde) momentum there. The Rossby wave's influence 
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Figure 11. Same as Fig. 10 at the 700 hPa level. 
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Figure 12. Eddy momentum flux cospectra at 400 (left) and 700 (right) hPa in the statistically steady 
state. Note that the global wave now has an equatorial phase speed of ~5 m/s. 



on equatorial winds - and in particular its deceleration - 
is diminishing in this phase, since as equatorial winds 
accelerate the wave no longer experiences a low-latitude 



critical layer. By day 250, this wave ceases to produce a 
momentum flux convergence at the equator. 

In each of the four snapshots of Fig. 10, a faster global 
wave with an equatorial phase speed of 12 m/s domi- 
nates the cospectrum at low latitudes. At the equator, 
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Figure 13. Sequence of snapshots of zonal wavenumber one, —0.6 days -1 frequency geopotential (shaded) 
and potential vorticity (contours) with zonal mean zonal wind multiplied by a factor of 5 overplotted 
(line) at the 400 hPa level for the Rot = 10.5 case. During the course of the transition to superrotation, 
particularly days 210-310, note how the disturbance remains coherent in latitude even as it propagates 
eastward (as is evident in Fig. 10) and that it leans into the horizontal shear. Note that the wave travels 
once around the globe in 1/0.6 ~1.7 days, and therefore eastward propagation of the pattern is not evident 
due to aliasing. 



this mode travels eastward relative to the mean winds 
from days 230 to 270. The momentum flux of the global 
wave is convergent at the equator, depositing eastward 
(prograde) momentum there. The wave maintains coher- 
ence past what appears to be a critical layer at 40° N/S 
latitudes, and as a result, the global wave propagates 
westward relative to the mean flow poleward of these 
latitudes. The eddy momentum flux is oppositely ori- 
ented in the westward-propagating component so that 
the flux divergence decelerates the mean winds at 40° 
N/S latitude. Between days 250 and 290, the zonal 
mean zonal wind at the equator is more than doubled 
by the momentum convergence of the global wave. By 
day 290, mean equatorial and mid-latitude winds have 
accelerated beyond the phase speed of the equatorial 



portion of the wave and the mode responsible for equa- 
torial acceleration no longer converges momentum to 
the equator. At this point, the mean flow adjusts due 
to the transport of westward (retrograde) momentum to 
the high-latitude jets, as can be seen in the diminished 
winds at 70° N/S latitude. 

Fig. 11 displays the eddy momentum flux cospectra at 
the 700 hPa level for the same sequence of times during 
spinup shown in Fig. 10. The same global disturbances 
dominate the eddy momentum fluxes at this level, with 
some important differences. The faster wave converges 
momentum to the equator as it does at the 400 hPa 
level, but it does so globally without a reversal in the 
direction of momentum flux (as occurs at the 400 hPa 
level poleward of 40 N/S latitudes). Momentum fluxes 
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at the 700 hPa level due to the slower wave are also 
oppositely oriented relative to those at the 400 hPa level. 
The 700 hPa level resides at the top of our boundary 
layer, and a frictional torque prevents equatorial superro- 
tation from developing despite strong eddy momentum 
flux convergence there. 

Fig. 12 displays the 400 and 700 hPa (left and right, re- 
spectively) eddy momentum flux cospectra in the statisti- 
cally steady state, i.e., after superrotation is established. 
A single, global wave dominates the cospectrum at a 
slower phase speed, ^5 m/s at the equator, than occurs 
during active spinup. The global wave at the 700 hPa 
level deposits prograde momentum at the equator. The 
equatorial eddy momentum flux is not well-organized 
at the 400 hPa level during this phase. However, the 
free atmosphere at this level does not directly experi- 
ence frictional torques from boundary layer processes 
(which are confined to p > 700 hPa) , has no significant 
zonal mean overturning circulation (as can be seen in 
the upper left panel of Fig. 2), and does not experience 
a vertical divergence of eddy activity (as can be seen in 
the left panel of Fig. 7). Consequently, the relatively 
weak and disorganized horizontal eddy momentum flux 
convergence is able to maintain the superrotating state. 
Section 5.5 contains further discussion of the mechanism 
maintaining superrotation. 

5.3. Structure of the global wave producing 
superrotation 

The spectral diagnostics in Figures 9-12 together 
reveal the disturbance responsible for equatorial momen- 
tum convergence during spinup has zonal wavenumber 
one and frequency -0.6 days -1 . Fig. 10 also shows the 
disturbance propagates both eastward and westward 
relative to the zonal mean zonal wind, depending on 
the latitude, but, importantly, the disturbance retains 
a coherent latitudinal structure. Figure 13 displays 
a sequence of snapshots of the 400 hPa geopotential 
(shaded) and potential vorticity (PV, contours) during 
spinup filtered for zonal wavenumber one and frequency 
—0.6 days -1 (the zonal mean zonal wind multiplied by 
a factor of 5 is also displayed). 2 During days 130 to 
210, the equatorial (eastward-propagating) component 
is beginning to phase with the high-latitude (westward- 
propagating) component. Both components strengthen 
during days 230 to 290, during which time the equato- 
rial wind speeds double in magnitude due to horizontal 
eddy momentum convergence there. The components 
then weaken to roughly their strength at day 130 and 
de-phase, thus ending a period of rapid acceleration of 
equatorial zonal winds. 

Near the equator, the global wave at times seems 
to resemble an equatorial Kelvin (gravity) wave with a 
characteristic chevron shape. However, the phase speed 
of the equatorial portion (12 m/s) does not correspond 
to typical gravity wave speeds, the gravest of which 
would have c >100 m/s. Poleward of 40 N/S latitudes 



(the location of the critical layer in Figure 10) the wave 
propagates westward relative to the mean flow but it 
remains anti-phased with the equatorial portion as it 
evolves. 

Fig. 14 displays the geopotential height anomaly asso- 
ciated with wavenumber 1 disturbances contributing to 
eddy momentum flux on day 270 (during spinup) with 
a frequency of -0.6 days -1 . Contours are spaced evenly 
from -0.5 to 0.5 m in increments of 0.1 m. The top row of 
panels show the pressure-longitude cross-section at three 
latitudes (the pattern moves toward the right in these 
figures). The disturbance is evidently the first baroclinic 
mode, as can be most clearly seen in the vertical slice at 
80N latitude (top-right panel). At lower latitudes, the 
vertical structure is more complicated which may be the 
result of shearing by the zonal mean zonal wind. The 
bottom row displays latitude-longitude cross-sections of 
the wave at the pressure levels indicated (for reference, 
the lower left panel corresponds to the level displayed in 
Fig. 13). At 900 hPa, the wave appears to be coherent 
across all latitudes with a slight northwest-southeast 
(northeast-southwest) tilt in the northern (southern) 
hemisphere. This tilting is more pronounced at 700 hPa, 
and still more so at 400 hPa. 

Fig. 15 displays the geopotential height anomaly as- 
sociated with wavenumber 1 disturbances contributing 
to eddy momentum flux in the statistically steady-state 
(day 1430 is shown). This disturbance, with a frequency 
of -0.25 days -1 , is clearly barotropic aside from mod- 
erate tilting that is evident in the frictional boundary 
layer. Its horizontal structure is very uniform across 
all latitudes except for in the frictional boundary layer 
where the contours tilt northwest-southeast (northeast- 
southwest), in the same manner they did for the wave 
contributing to spinup in Fig. 14. 

The cause of the northwest-southeast (northeast- 
southwest) tilting and its influence on the zonal mean 
zonal wind is now discussed. 

5.4. Structure of the zonal mean potential vorticity 

The development of the instability can be seen in the 
time-evolution of zonal mean PV and zonal winds in 
Fig. 16, shown at the 400 hPa level of the Ro T =10.5 
case on the days indicated. A necessary condition for 
bartotropic instability is that the meridional gradient of 
zonal mean PV, dq/dy, change sign in the horizontal. 
At day 180 of Fig. 16, mean equatorial winds are very 
calm and the equatorial PV is nearly zero aside from 
small perturbations. Zero PV is an indication of an 
angular-momentum-conserving overturning circulation 
during the initial, nearly axisymmetric phase of the sim- 
ulation. The PV gradient at the equator has steepened 
by day 200, which introduces areas of negative PV gradi- 
ents near 20° N/S latitude. Equatorial winds accelerate 
during this time, which is suggestive of a link between a 
large-scale barotropic instability and the emergence of 
superrotation. 
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Figure 14. Vertical and horizontal slices of zonally asymmetric geopotential filtered for zonal wavenum- 
ber 1 and frequency —0.6 days - (shaded) during spinup (day 270) at the latitudes or levels indicated. 
Contours are spaced evenly from —0.5 to 0.5 m in increments of 0.1 m. 
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Figure 15. Same as in Fig. 14 for the global wave, now with frequency 
steady-state (day 1430). 



-0.25 days 1 , in the statistically 



The zonal mean PV gradient in the Rot = 10.5 case 
during spinup is displayed in left panel of Fig. 17. The 
zero PV gradient line is marked with a magenta contour, 
and three isotherms are over-plotted in black. There is 



clear evidence for a PV gradient reversal near the surface 
and at high latitudes, indicating a role for baroclinic 
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Figure 16. Same as in Fig. 4 for snapshots during spinup (on days indicated) of the Rot = 10.5 case at 400 hPa. 
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Figure 17. (Top row) Zonal- and time-mean PV gradient, dq/dy, during an active phase of spinup (left 
panel) and in the steady-state (right panel) for the Rot = 10.5 case. The zero PV gradient line is marked 
with a magenta line, with warm (cool) colors indicating positive (negative) PV gradient. The 295K, 300K, 
and 305K isentropes are overplotted in black. (Bottom Row) Same on isentropes. White lines indicate 
areas of large, positive PV gradient. 



instability during spinup. The global (first) baroclinic 
mode responsible for spinup (Fig. 14) is associated with 
this high-latitude baroclinic instability. Being global in 
scale, the disturbance is unable to propagate either ver- 



tically or horizontally and thus our notions based on the 
propagation and breaking of linear waves is insufficient 
to explain its interaction with the mean flow. 

The PV gradient also reverses sign on either side of 
the equator at the 400 hPa level and at 30 N/S latitudes. 
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Figure 18. Schematic of barotropic instability during 
spimip of the Rot = 10.5 case. The meridional structure 
of the zonal mean zonal wind is shown as the line. In- 
teracting edge waves form circulation patterns that slant 
into the shear of the zonal winds. The resulting Reynolds 
stresses converge momentum to the equator. 

Eastward-propagating Rossby waves initiated in these 
regions of PV reversal span the equator and interact 
with one another. This accounts for the slow, eastward 
propagation (relative to the gravity wave phase speed) of 
the equatorial portion of the global wave. The eastward 
propagation of the equatorial portion together with an 
eastward, horizontal zonal wind shear in zonal winds al- 
low phase-locking between high-latitude (westward) and 
low-latitude (eastward) portions of the global wave (as 
is a characteristic of barotropic instability - see Section 
5.5). 

The zonal mean PV gradient in the Rot = 10.5 case 
in the statistically steady state is displayed in right panel 
of Fig. 17. The regions of negative PV gradients at the 
high-latitude surface are no longer present, indicating 
an absence of baroclinic instability. The high-latitude 
PV gradient reverses sign in the horizontal, indicating 
barotropic instability is associated with the barotropic 
mode seen during the steady-state (Fig. 15). The middle 
troposphere at low latitudes no longer has PV gradient 
reversals that were present during the spinup phase, al- 
though there are small regions of negative PV gradient 
confined to the boundary layer (pressures greater than 
700 hPa). 

The PV diagnostics in Figs. 16 & 17 together sug- 
gest a role for baroclinc and barotropic instabilities in 
the emergence and maintenance of superrotation in our 



simulations. It is unclear from our experiments whether 
both instabilities are necessary for the transition to oc- 
cur. The mechanism giving rise to a convergence of 
eddy momentum flux at the equator is very similar in 
character to barotropic instability, which we now discuss 
in more detail. 

5.5. Barotropic instability 

The time-evolution of geopotential and potential vor- 
ticity shown in Fig. 13 and Fig. 16 and the reversals 
of PV gradients present in Fig. 17 are both suggestive 
of a role for barotropic instability during spinup of the 
superrotating state. Fig. 18 displays a schematic of 
the mechanism of barotropic instability for an ideal- 
ized zonal mean zonal wind following the initial, nearly 
axisymmetric phase of the simulation. Rossby waves 
form in regions of vorticity gradients and if the vorticity 
gradient is strong these are essentially edge waves. If 
a pair of edge waves form at the equator and at higher 
latitudes and have sufficient overlap in space, they are 
able to phase lock and interact to form an instability. 
This gives rise to the circulation depicted in Fig. 18 by 
slanted ellipses, with the instability 'leaning into the 
shear' in order to extract energy from it. Since the zonal 
wind speed is higher at high latitudes than at low lati- 
tudes, phase locking of the waves requires them to travel 
in opposite directions relative to the local mean flow, 
and for Rossby waves this would require their local PV 
gradients (dq/dy) to be of opposite sign. Because the 
circulation formed by the instability slants into the shear 
the Reynolds stresses from these disturbances transfer 
eastward momentum down the wind shear, converging 
momentum at the equator and producing superrota- 
tion. Furthermore, because the interacting edge waves 
must span a large meridional distance, the zonal scale 
of the waves is correspondingly large, so explaining the 
emergence of wavenumber one in our simulations. 

The above picture is, however, an oversimplification, 
even if it captures the basic underlying mechanism lead- 
ing to momentum convergence at the equator. At the 
onset of the instability (for example, at day 220 in Fig. 
16), the potential vorticity gradient is positive at the 
equator and becoming negative in the subtropics before 
changing sign again to be large and positive at high 
latitudes. These large scale reversals in potential vortic- 
ity gradient with latitude are certainly indicative of a 
barotropic instability, although it is not only the hori- 
zontal shear and the beta term that contribute to the 
potential vorticity; rather, terms involving the strati- 
fication and vertical shear are equally important (not 
shown). 

Given the large-scale nature of the instability, the 
reason for the disruption of the development of superro- 
tation by a sponge layer placed at a particular latitude 
band becomes clear (see Fig. 6). The structure of the 
global disturbance is significantly altered if any portion 
of it is damped, thus compromising the wave's ability to 
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Figure 19. Same as Fig. 10 for the Rot = 1.3 case during spinup. 



create momentum flux divergence. We found the disrup- 
tion to be particularly effective for a thin sponge layer 
placed at the equator, consistent with an equatorial role 
in the instability giving rise to superrotation. 

We can also see what prevents strong superrotation 
from developing in the Rot = 1.3 case. Fig. 19 displays 
a sequence of phase-speed-latitude cospectra (shaded) 
with the zonal mean zonal wind overplotted (line) for 
the Rot = 1.3 case. At high latitudes, several Rossby 
waves that propagate westward relative to the mean flow 
are evident. These disturbances are due to baroclinic 
instability, and they converge westward (retrograde) an- 
gular momentum to critical layers throughout the entire 
latitudinal domain. A global wave transiently appears 
and accelerates the equatorial winds for a time, but the 
restoring torque of the high-latitude waves breaking at 
lower latitudes dominates the low-latitude momentum 
budget, keeping equatorial winds weakly prograde. 

6. Superrotation in the Solar System 

The atmospheres of two of the four terrestrial bodies 
in the Solar System with thick atmospheres — Titan 
and Venus — are observed to superrotate, whereas the 



atmospheres of Earth and Mars do not. We have shown 
the thermal Rossby number is a control parameter in the 
transition to superrotation, and this parameter is most 
sensitive to rotation rate and radius. Titan and Venus 
respectively have rotation periods 16 and 243 times that 
of Earth; Venus is nearly Earth's twin in size, while 
Titan's radius is ^40% Earth's, and consequently the 
thermal Rossby numbers of Titan and Venus are both 
much greater than one. Since the transition to superro- 
tation occurs for values of Rot exceeding one, Titan and 
Venus lie well beyond this transition in parameter space. 
It is possible that Venus (with Rot ~ 1200) is in yet an- 
other parameter regime from that explored in this paper, 
but the essential dynamics of the superrotation of Titan 
(with Rot 30) may be similar to those occurring in our 
simulations (this connection, and the role of seasonality, 
will be further explored in a subsequent paper). The ra- 
dius of Mars is roughly 50% of Earths and has a nearly 
identical rotation rate. Therefore, iZox(Mars) ~ 0.1, 
and so Mars is the closest to being at the superrotation 
transition; a further reduction in the Martian rotation 
rate could cause the atmosphere to superrotate. How- 
ever given Mars's lack of a large satellite and weak tidal 
coupling with the Sun, it is unlikely that the rotation 
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rate has changed significantly over the lifetime of the 
Solar System (or will change in the foreseeable future). 
Supcrrotation on the gas giants seems likely via mecha- 
nisms that do not involve a global mode but that may 
involve convection at the equator [Yamazaki et at, 2005; 
Lian and Showman, 2009; Schneider and Liu, 2009] 
or deep convection [Heimpel and Aurnou, 2007; Kaspi 
et al, 2009]. 

7. Conclusions 

We have described the mechanism of transition of 
a terrestrial atmosphere to superrotation that occurs 
when Rot, the thermal Rossby number, exceeds unity, 
as is characteristic of the atmospheres of Venus and Ti- 
tan. The Earth's equatorial atmosphere (with Rot <C 1) 
does not superrotate because there is no mechanism 
to excite the large-scale disturbance we have identified 
in our Rot = 10.5 case. Rather, Rossby waves that 
are small compared to the domain size are formed at 
mid-latitudes, propagate equatorward and break, decel- 
erating the equatorial flow and damping any tendency 
toward superrotation that might be present. However, 
if Rot exceeds unity a robust mechanism emerges: the 
Rossby wave source in midlatitudes weakens (because 
the scale of baroclinic instability becomes larger than 
the planetary circumference) and a global disturbance 
predominantly of zonal wavenumber one develops that 
provides a convergence of momentum at the equator. 
Both baroclinic and barotropic instability play a role in 
the spinup phase, while barotropic instability alone is 
responsible for the maintenance of superrotation in the 
statistically steady-state. 

The global disturbance during spinup of our Rot = 
10.5 case as the first baroclinic mode is associated with 
high-latitude PV gradient reversal in the vertical di- 
rection (suggestive of baroclinic instability) as well has 
PV gradient reversals in the horizontal direction, near 
the equator (suggestive of barotropic instability). The 
disturbance deposits eastward (prograde) momentum 
in the equatorial region, accelerating the zonal wind 
into a superrotating state. The PV gradient changes 
sign in the horizontal direction near the equator, indi- 
cating a role for barotropic instability. Disturbances 
associated with these near-equatorial regions overlap, 
giving the equatorial portion of the mode the appear- 
ance of a chevron that is symmetric about the equator. 
The spinup mechanism thus involves a mixed baroclinic- 
barotropic instability; the relative importance of these 
needs further study. 

Superrotation is maintained in the steady-state by 
a persistent, global barotropic mode associated with a 
mixture of high- and low-latitude regions of PV gradient 
reversal in the horizontal (indicating barotropic insta- 
bility). Equatorial eddy momentum flux convergence 
by the global barotropic mode is generally disorganized 
and weak. However, frictional and other torques tend- 



ing to decelerate mean zonal winds at the equator and 
above the boundary layer are also quite weak, so little is 
required of the global disturbance in order to maintain 
the superrotating state. 

The heuristic picture of linear wave propagation and 
critical layer absorption is in contrast with the global 
coherence of the dominant waveforms in our high-Rossby- 
number simulations. The disturbances responsible for 
the transition to superrotation are global in scale and 
thus are unable to propagate. Rather than providing 
a critical layer for a propagating wave, the zonal mean 
latitudinal shear plays an important role in the merid- 
ional transport of momentum by its interaction with the 
global wave. 

The dominance of long-wavelength disturbances in 
our results suggest modeling studies of Titan aimed at 
developing realistic levels of superrotation should not 
require high resolution. On the other hand, we demon- 
strated the development of superrotation in our model 
is strongly sensitive to damping at the equator. Models 
having unrealistic dissipation at or near the equator will 
thus have difficulty producing and maintaining superro- 
tation. 

For intermediate Rot, there is a competition between 
equatorial acceleration by global waves and decelera- 
tion by Rossby waves due to mid- and high-latitude 
baroclinic instability. Our Rot = 1.3 case shows evi- 
dence for transient global waves like that seen in the 
Rot = 10.5 case, but also displays baroclinic instability 
at mid-latitudes. Strong superrotation does not develop 
because of the deceleration caused by the breaking of 
meridionally propagating Rossby waves generated at 
mid- and high-latitudes by baroclinic instability. 

Equatorial waves with a characteristic chevron shape 
similar to the equatorial portion of the global wave in our 
simulations are observed in Venus' cloud layer, and these 
waves may be connected with the mechanism described 
here. If the observed waves correspond to the equatorial 
component of the global wave we have described, then 
a high-latitude component should also be present whose 
geopotential anomaly varies coherently with respect to 
the equatorial component (see Fig. 13). Similar waves 
have not yet been observed in Titan's atmosphere, but 
may nonetheless be observable from time to time in 
thermal infrared data from the Cassini orbiter. 

Of course, notwithstanding the heuristic description 
that we have given, the mechanism we have identified 
remains to be analytically described in terms of a classi- 
cal linear instability or wave analysis that would yield 
its spatial structure and dispersion relation in terms of 
the fundamental parameters of the flow, in a manner 
akin to the classical analyses of the eponymous waves 
and instabilities of Rossby, Charney and Eady. 
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Notes 

1. Throughout this paper, we use the word superrotation to de- 
scribe an atmosphere that has more axial angular momentum 
than the solid planet at the equator. This will usually occur 
as an angular momentum maximum at the equator, because if 
the atmosphere is supcrrotating at higher latitudes and not at 
the equator it will be inertially unstable. This meaning is to 
be contrasted with 'global superrotation', which occurs when 
the integrated angular momentum of the atmosphere is more 
than it would be in solid body rotation. Global superrotation 
is not uncommon in planetary atmospheres, and occurs in the 
Earths atmosphere. 

2. Filtering is performed by direct FFT in longitude and time, se- 
lecting the desired wavcnumbcr and frequency, and performing 
a reverse FFT on the filtered spectrum. 
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